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Abstract 

We consider spectral semi-Galerkin approximations for the strong so- 
lutions of the nonhomogeneous Navier-Stokes equations. We derive an 
optimal uniform in time error bound in the H 1 norm for approximations 
of the velocity. We also derive an error estimate for approximations of the 
density in some spaces L r . 

1 Introduction 

Let il C R™, n = 2 or 3, a C 1 ' 1 — regular bounded domain, and T > 0. We are 
interested in the initial boundary value problem 

' p(u t + u ■ Vu) - Au + Vp = pi in fix [0, T), 
divu = in fix (0,T), 

Pt + u-Vp = in fix (0,T), n n 

' u = on dQ x [0,T), 1 ' 

u(x, 0) = u Q (x), e fi, 
p(x,0) = p (x), \fx e tt. 

These are the equations of motion for nonhomogeneous incompressible fluids. 
The unknowns are the velocity u(x, t) £ E™ of the fluid, its density p(x, t) 6 R 
and the hydrostatic pressure p(x,t) G R. The functions uo(x) and po(^) are 
respectively the initial velocity and initial density. The function i(x, t) is the 
density by unit of mass of the external force acting on the fluid. Here, without 
loss of generality to our aim, the viscosity is considered to be one. 
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Different methods have been used to study the existence, uniqueness and 
regularity of solutions for problem ifTTjlfsee g|, [3], 0, [SI 0, QUI, EH, G2, 
H3|, U, E3, EH)- Here, we are interested in the spectral Semi-Galerkin 
method applied to system jl.ip. the word spectral being used to indicate that 
the eigenfunctions of the associated Stokes operator are used as a basis of ap- 
proximations. It is very important to derive error bounds for Galerin methods, 
due to the wide application of these methods in numerical experiments. Even 
the case of spectral Galcrkin method may be used as a preparation and guide 
for the more practical finite element Galerkin method. 

A systematic development of error bounds for the spectral Galerkin method 
applied to the classical Navier-Stokes equations was given in ^j- These error 
bounds are local in time in the sense that they depend on functions which grow 
exponentially with time. As observed in jS], this is the best one may expect 
without any assumptions about the stability of the solution being approximated. 
Optimal uniform in time error estimates for the velocity in the Dirichlct norm 
were also derived in [S], assuming uniform boundedness in time of the L 2 — norm 
of the gradient of the velocity and exponential stability in the Dirichlet norm of 
the solution. 

For the variable density case, error bounds were first obtained in |16| . where 
local in time error bounds were derived. Moreover, a result of uniform in time 
error estimates in the L norm was also stated. This last result, however, is not 
optimal. Moreover, it requires the assumption u £ L°°(0, T, H 3 (£!)). As pointed 
out in [7], this assumption is pretty restrictive, since it requires a global compati- 
bility condition on the initial data even for the classical Navier-Stokes equations. 
Error bounds were also derived in 0] without explicitly assuming stability, but 
requiring exponential decay of the external force field. This hypothesis though 
is very restrictive as well, since gravitational forces do not satisfy it. 

Here we derive error bounds assuming the solution (u, p) to be po- conditionally 
asymptotically stable, a notion to be defined in section [21 The number po is re- 
quired to satisfy 6 < p < oo, and is related with the regularity of allowed 
perturbations of the density equation. In a similar notion has been used to 
treat the classical Navier-Stokes equations(see also EH]). Here, we adapt it in 
the proper way to be used in the variable density case. With this assumption, 
we obtain an uniform in time optimal error bound in the Dirichlet norm for the 
velocity. An error bound depending on time for the density in some spaces L r 
is also derived. 

In sectionElwe state some preliminary useful results. In section[2|we describe 
the approximation method, the stability notion to be used, and state the main 
result. Section 0] is dedicated to deriving a priori estimates. Finally, we present 
the proof of the main result in section [5] 

To simplify the notation, we denote by C a generic finite positive constant 
depending only on fl and the other fixed parameters of the problem. As usual, 
it may have different values in different expressions. When necessary, we em- 
phasize that the constants may have different values using the notation C\, Ca, 
and so on. 
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2 Preliminaries 



Throughout this work, we consider the usual Sobolev spaces 

W m >«(D)={feL«(D),:\\d a f\\ L9{D) <+oc,]a|<m}, 

for a multi-index a, a nonnegative integer m and 1 < q < +00, where the 
domain is D = il or D = ft x (0, T), < T < +00, depending on the context. 
We write H m (D) := W m > 2 (£>) and denote by H™{D) the closure of C °°(L>) 
in H m (D). If P is a Banach space, we denote by L q ([0,T];B) the Banach 
space of P-valued functions defined on the interval [0, T] that are L 9 -integrable 
in Bochner's sense. Spaces of M™ valued functions, as well as their elements, are 
represented by bold face letters. We write 

C^ CT (ft) := {v e Cg°(ft)/ div v =0 in f2}, 

and denote by H and V the closure of C§° .(f2) in L 2 (il) and H 1 (fi) respectively. 

Throughout this work, the orthogonal projection from L 2 (f2) onto H is writ- 
ten as P. Thus, the well known Stokes operator is —PA, defined on VnH 2 (£l). 
Its eigenf unctions and eigenvalues are denoted by w fc and respectively. The 
usual L 2 (f2) inner product and norm are respectively indicated by (•, •), and || • ||. 

It is well known that {w fc (a:)}^ 1 form an complete orthogonal system 
in the spaces H, V and V n H 2 (f2) equipped with the usual inner products 
(u, v), (Vu, Vv) and (PAu, PAv) respectively. 

For each k e N, we denote by Pk the orthogonal projection from L 2 (S7) onto 
V fc := span[w\ w fc ]. For all f, g G L 2 (£l) and k,m <E N, it holds 

(i) (P fc f,g) = (f,P fe g), 

(*») (Pf,g) = (f,Pg), 
(Hi) ((P m - P fc )f, g) = (f, (P m - P fe )g), 

(iv) ((P-P fc )f,g) = (f,(P-P fc )g). 
The following bounds are useful to our ends. 
Lemma 2.1 (Rautmann[15j). IfvGV, then 

||t>-J\«|| 2 < J-||V»|| a . 

Moreover, if v e V n H 2 (S1), i/ien 

||«-PH| 2 <t^||PAu|| 2 , 
A fe+1 

||V«- VJ-fc^l 2 < -l-||PAt;|| 2 . 
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Remark 2.2. From Lemma it follows that if f G H 1 (ri), then 
IKl-^Pfll^y^llVPfjl 2 . 

Moreover, since P : H 1 (J7) — > H 1 (J7) is a continuous operator^], we have 

||VPf|| 2 <C||f||^. 
Thus, for all f G H^fi), one has 

\\{I-Pk)Pf\\<^-\\^ 

Since PPk = PkP = -Pfe, one equivalently obtains 

HPf-Pfeff^^llfll^. 

Moreover, the above relations also hold if one replaces P by any P m , m > k. 
Analogously, one may check that 



||(/-P fe )Pf|| 2 <T^||f||k 
A fc+1 



for all f G H 2 (f2). We also have, as an easy consequence of the L 2 -orthogonality 
of the functions {w fe }^° =1 , the following: Let m > k, m, k £ N, f G L 2 (fi) and 
v m G V m , v fe G V fc . Then 

((P m - P k )f, v m - v fc ) = (f, (/ - P fc )v m ). 

3 Stability concept and main result 

We consider problem HI. 1(1 f° r a h time t > 0, and suppose the data given therein 
satisfy 

u GVnH 2 (f7), (3.1) 
sup ||f|| H i < oo ; sup ||f t || < oo, (3.2) 

t>0 t>0 

po G C\Q) ■ < a < Po < /3, (3.3) 

where a and (3 are constants. We also suppose that there exists M > such 
that the solution (u, p) of 1(1.1(1 satisfies 

||Vu(t)|| < M ; Vt > 0. (3.4) 

If ra = 2, then conditions ((3.1(1 and ((3.21) imply that 1(3.4(1 holds. If n = 3, then 
inequality 1(3.4(1 holds for f and Uo sufficiently small (see |S]). For a given po 5 
6 < Po < °°7 we also assume (u, p) to be po- conditionally asymptotically stable 
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(see jHlEl f° r similar notions of stability). To define this notion of stability, we 
first define perturbations of system Ql.l|l . The pair (£(x,t),r)(x, i)), defined for 
t > to > 0, is called a perturbation of (u, p) at time to if (u := u+£, p := p+rj) is 
a solution of i|l.lfl . within = 0. Therefore, setting £ :— £(-,t ), i] := 77 ( - , to ) , 
the pair (£, p) is a solution of the initial boundary value problem 

' n t + ?(u • v)$ + m ■ v)u + p(t • v)« + v q = 

= A£ + (p - p)(u t + (u • V)u - f) in n x (t , 00), 
ft + ((£ + u)-V)p = in fix (t , 00), 
< V ■ £ = in fi x (t , 00), 
Z\dn = 0, 

£(s,*o) = ^oW) Vx e ^> 

p(x,t ) = p(x,to) +r/ Q (x), Vx e n. 

,(3.5) 

Now, for a given po, 6 < po < oo 5 we define the concept of po- con( htional 
asymptotic stability. 

Definition 3.1. The pair (u, p) is said to be po-conditionally asymptotically 
stable if for all to > there exist positive numbers A, B, 5, M\, M2 and a 
continuous decreasing function F : [0, 00) — > ]R + , F(0) = 1, lim F(t) — such 

t — >oc 

that, for all £ G H 2 (n), m £ L°°(fi) n W 1,Po (fJ), satisfying ||V£ || < 5, 
||PA£ || < A, I] 770 1| < B, problem 1)3. 5fl is uniquely solvable with 

$ eLL([to^);VnH 2 (0)), 
^eLf^Ito.ooJjH 1 ^)), 

77 e L°°([to, oo);i°°(n) n w^^fi)). 

Moreover, 

||V?7(-,t)|| LP o < Mi , V* > to, (3.6) 
||V£(-,t)|| <M 2 ||V^ ||F(t-i ), Vt>t . (3.7) 

Remark 3.2. We use a general function F(t) in Definition 13. II just to stress out 
that the results here do not require an exponential decay rate. 

The solution of problem can be obtained through a spectral semi- 

Galerkin approximation, that is, a spectral Galerkin approximation 

n 

U n (x,t)=Y,Ckn(t)^ k (x) 
k=l 

for the velocity u, uniquely determined by 

(p n ul\ cp n ) + (p n u n ■ Vu», 0") + (Vu™, cp n ) = (p n i, cp n ), t > 0, (3.8) 
(u"M)- 110(00,0") = 0, (3.9) 
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for all <p n of the form 4> n (x) = ^""^ dk"w k (x), and an infinite dimensional ap- 

fc=i 

proximation p n for the density, solution of 

P n (0) = po. 

It can be proved that (u", p n ) converges in an appropriate sense to (u, p), solu- 
tion of Our main result is 

Theorem 3.3. Suppose {u, p) to be po- conditionally asymptotically stable, for 
some po, 6 < po < 00 • Then, there exist constants N 6 N and C > such that 
if n> N then, for all t > 0, 

||V«(-,t)-Vu"(.,t)||<— - ^-r. (3.11) 

Moreover, if 6 < po < oo, i/ien 

||p(.,t)-p»(.,f)||^<^t_ 2 <r<-^, (3.12) 

and if po = 00, i/ien 

Ct 

\\p(;t)-p n (;t)\\ L r<- -r , 2<r<6. (3.13) 

(A„+l)2 

TTie constants N , C, depend only on the domain, on the norms of the data in 
fff.ij) . Q3. 2\) and on the constants introduced in \3.4\j and definition \3.1\ 

4 A priori estimates 

We first state a general simple result that will be used later on. A proof is given 
in the Appendix. 

Lemma 4.1. Let h(t) be an integrable nonnegative function. Suppose there 
exist nonnegative constants a\, 0,2 satisfying 

t 

h(r)dr < d\(t — t ) + a 2 , 



for all t, to with < to < t. Then, 

sup I e~* / e T h{r)dT I < 00. 
t>o \ Jo / 

For u, solution of 1)1. If) , and the perturbations we have: 
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Lemma 4.2. Given e > 0, the bounds 

i|||V M (-,t)|| 2 + ^||PA<,i)|| 2 + i|| p ^(.,t)|| 2 < (4.1) 

C(||/M)|| 2 + ||Vu(.,t)|| 6 ), 

^ t( '' t)l12 + ^ l|V ^'' t)l12 + ^H PA ^-' f )H 2 ^ ^ 
C(\\Vt(;t)\\ 2 \\PAu(;t)\\ 2 + ||V£(-,i)|| 6 + K(-,t)|| 2 + \\PAu(-,t)\\ 2 ) , 

hold for all t > 0. 

Proof. Inequality 1)4. 1|) was proved in ^Sj. Inequality l|4.2[) can be proved in a 
completely analogous way. □ 

Corollary 4.3. For all t>to, we have 



f \\PAu(-,T)\\ 2 dT<C + C(t-t ), 
Jt„ 



(4.3) 



one gets 



/ \\ut(;T)fdT<C + C(t-t ), 
J to 


(4.4) 


f \\PAt(:T)\\ 2 dT<C+C(t-t ), 

J to 


(4.5) 


[ M t (-,T)\\ 2 dT<c + c(t-t ). 

J to 


(4.6) 


combining inequalities \4-Sfy, \4-4\li \4-®il an ^ \4-ty 


wrei/i lemma\4~^\ 


supe"*/ e T \\ut(-,T) \\ 2 dr < oo, 
t>o Jo 


(4.7) 


supe"*/ e T \\ P Au(-, t)\\ 2 dr <oo, 
f>0 Jo 


(4.8) 


supe"* / e T £ t (-, r) 2 dr < oo, 
t>0 Jo 


(4.9) 


supe"*/ e T ||PA|(-,T)|| 2 dr < oo. 
t>o Jo 


(4-10) 



The following lemma states some bounds for u which are very important for 
our later arguments. 

Lemma 4.4. We have 

sup||«t(-,t)|| < oo, (4-11) 
t>o 

sup||PAu(-,t)|| < oo, (4.12) 

t>Q 

supe"'/ e T ||V%(-,r)|| 2 dr < oo. (4.13) 
t>o Jo 
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Proof. We begin by proving (|4.12[) . supposing l|4.11|l holds. Setting v = —PAu 
in the weak formulation of problem , we get 



Thus, 



- (pu t , PAu) - (pu ■ Vu, PAu) + 1 1 PAu| | 2 = - (pf, PAu) . 



|PAu|| < ||pu t || + ||pu.Vu|| + ||pfl| 

< /3||u t ||+/?||u|| LS ||Vu|| L3 +/3||f|| 

< /3||u t ||+C^||PAu||5||Vu||i+/3||f|| 



< /3(l|ut|| + ||fj| 



^||Vu|| 3 + i||PAu|| 
1, 



< 2/?||u t || + 2p\\i\\ + (CpYM* + -||PAu||. 

Therefore, by l|4.11|l and l|3.2[) . we have 

sup||PAu(-,t)|| < 2/3 sup ||u t (-,t)|| + 2/3 sup ||f(-,i)|j + (C(3) 2 M 3 < oo, 

t>0 t>Q t>0 

which proves (|4.12(l . To prove (|4. 1 1|) and (|4. . differentiate the weak formu- 
lation of problem , and set v = u t to get 

(p t u f , ut) + (pu u , u t ) + ((p(u • V)u) t , u t ) + (Vu t , Vu t ) = ((pf) t , u t ). 

Therefore, 

^ll^ut|| 2 + l|Vu 4 || 2 = ~J P U-V(VLfU t )dx 

-((p(u-V)u) i)Ui ) + ((pf) t) u t ). (4.14) 
Now, bound each term on the right hand side of (|4.14() as follows: 



~ / pu- V(u t • u t )dx 
z Jo. 



< ||/0||o C ||u|| L 4||u t || L 4||Vu t || 

< C/3||Vu||||Vu 4 ||{||u t ||^||Vu t ||f} 
= C/3||Vu||||u t ||i||Vu t ||J 

< C e (C/3||Vu||) 8 |H| 2 + e ||Vu t || 2 , 



l(pf*,u t )| < 



if 



\m 2 



l(Ptf,u t ) 



p t f • u t da; 



/ div(pu)f • u t dx = — I pu ■ V(f • u t ) dx 

Jn Jq 



< ||p|U~||u||||V(f-u t )|| <C £ (C/3||Vu||||Vfl|) 2 + e||Vu t || 2 

< C e (C/3M||Vf||) 2 + e||Vu t || 2 , 



|(p(u,.V)u,u t )| < ||p||^||u t ||| 4 ||Vu|| <C/3||Vu||||u t ||T|| V u t p 

< ^(C/JllVullHuJ^ + ellVuJ 2 <C £ C 4 ^ 4 M 4 ||u t || 2 + e||Vu t || 2 , 

|(p(u-V)u t ,u t )| < ||p|Uoo||u t H * || Vu|| HVtxtH 5 < C* £ (C/3M) 8 ||u t || 2 + e||Vu t || 2 , 

\(p t (u- V)u,ut)| = | / pt(u- V)u-u t | = |- / div(pu)(u- V)u-u t | 
Jn Jn 

< CC e (3 2 M 4 ||PAu|j 2 +4e||Vu 4 || 2 . 
Therefore, from equation 14.14JI one obtains 

~|| VpM 2 + l|Vu t || 2 < C\(C(3Mf\\M 2 + yllfJ 2 + C e (C(3M || Vf||) 2 
+ C t C A (3 A M A \\u t \\ 2 + CC t f3 2 M 4 \\P Au|| 2 + ^8e + ||V Ui || 2 . 

^ ' (4.15) 

Now, choose e < — to conclude 
16 

lllv^u.^oip + qivu^^oi^c + qiPAui^^f + qiu^.,^!! 2 , 

(4.16) 

where C > is an absolute constant, and the constant C depends only on fi, 

||/9||l°o, ||Vf||, ||f t ||, sup||Vu||. Now, multiplying inequality (|4.16() by e* and 

t>o 

integrating over [0, i], one gets 

e*||Vpu t (-,t)|| 2 + C f e T \\Vu t (;r)\\ 2 dr < \\(y/pu t )(; 0)|| 2 + C [ e/\\PAu(-, r)\\ 2 dr 
Jo Jo 

+ C [ e T dr + C f e r ||u t (-,T)|| 2 dr. 
Jo Jo 

Hence, 

||>u 4 (-,i)|| 2 + C e ^ f e T \\Wu t (;T)\\ 2 dT < e-*/3||u t (.,0)|| 2 + Ce-* f e T \\PAu(-, r)\\ 2 dr 
Jo Jo 

+ Ce"* f e T dT + Ce~ t I e T \\u t {- , t)\\ 2 dr . 
Jo Jo 

Using inequalities (|4.7|l and (|4.8II . we get the desired result. 

Corollary 4.5. For all to, t G M, < to < one Zias 

rt 

\\Vut(-,T)\\ 2 dT < C(t-t )+C. (4.17) 



to 
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Proof. Integrating inequality i|4.16[l from to to t, one gets 

ft r t 



\\^pu t (;t)\\ 2 + C f ||Vu t (.,r)|| 2 dT < ||VM(.,t )|| 2 + C f dr 

J to J t 

C f \\PAu(;T)\\ 2 dr + C f ||u t (-,r)|| 2 dr 

J to J to 



< \\Jpu t (-,to)\\ 2 + C(t-t a ) 

+ C(t-t ) fsup||PAu(-,i)|| 2 



+ C(t-t ) sup||u t (-,*)|| 



t>0 



Using inequalities (|4.11|) and (|4.12|) . one obtains the bound l|4.17|l . □ 



A priori estimates for the solution £ of problem l|3.5|l . similar to those in 
Lemma 14.41 for u, also hold. Indeed, if ||V£ || < S, where 5 is the number 
referred to in Definition EH then it follows by IpTTjl that ||V£(-,*)|| < SM 2 . 
Therefore, the function u = u + £ is a solution of the nonhomogeneous Navier- 
Stokes equations satisfying ||Vu|| < M + 5M%. Moreover, if ||PA|(-, i )|| is 
bounded then ||PAu(-, to)\\ is also bounded. In this case, analogously to the 
proof of Lemma f4. 41 one can bound ||PAu(-, t)\\, for t > to. This bound implies 
that ||PA£(-,£)|| is bounded, for t > to- Summarizing, 

Lemma 4.6. For perturbations £ satisfying ||V£ || < 5 and ||PA£ | < C'o, 
C'o > 0, we have ||PA£(-,f)|| < C, for all t > to. The constant C depends 
on ||PA£ ||. C'o, fl, the initial data of problem and on the norms and 

constants appearing in \3. 61) and \3. 7| ). 

It also holds 

Lemma 4.7. The function £ satisfies 

f \\Vt t (;T)\\ 2 dT<C(t~to) + C, (4.18) 
Jto 

for all t , t, < t < t. 
Proof. Note that 

f \\VU,r)\\ 2 dr <c(f ||VS t (,r)|| 2 rfr + f || Vu t (, r)\\ 2 dr 

Jto \Jt Jt 

The second term on the right hand side of this inequality is bounded (Corollary 

I4.5fl . Therefore, it only remains to bound / |j Vut(-, T)\\ 2 dr. This bound follows 

Jt 

analogously to the bound for u. □ 
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Now, let u = Ck{t)w k (x) be the expression of u, the solution of Ijl.lfl . in 



fc=i 



terms of the eigenf unctions of the Stokes operator. Let v™ := Ck(t)w k (x) 



fc=i 



be its n-th partial sum, and define e" := u v™ and xp n :— u" — v™. We begin 
by bounding e™. 



Lemma 4.8. The bounds 

\\X7f> n (--f)\\ 2 < 

Ari+l 



Ve"(.,<)|| 2 < -°- (4.19) 



and 



||e"M)l| 2 < ^3— (4.20) 

hold for all t > 0. 
Proof. We bound 

||Ve"(-,t)|| 2 = ||V C k (t)w k (.)\\ 2 <—\\PA £ C k (t)w k (-)f 

k=n+l n+1 k=n+l 

< ^||PAu(.,t)|| 2 <-^, (4.21) 

^n+l -Ari+1 

as desired. Moreover, 

oo 1 oo 

||e»M)|| 2 = || J2 C k (t)w k (-)r<—\\X7 J2 C fc (t)w fc (.)|| 2 

k=n+l n+1 k=n+l 

< J-||Ve"(.,i)|| 2 <-^. (4.22) 

A n+l 

□ 

Now, we prove that a suitable bound for || X7xp n (-, t)\\ implies in a bound for 
||PA^(-,t)||. 

K 

Lemma 4.9. If for some constant K > the inequality || Vi/' n (-, t)\\ < 



holds on an interval [0,t*], then there exists C > 0, independent of n, such that 

\\PAip n (-,t)\\ 2 <C, (4.23) 

for all t G [0,t*]. 
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Proof: Since ip n = u™ v™ and sup ||PAv n (-, < sup ||PAu(-, t)\\ < oo, 

t>0 t>0 

one only needs to bound ||PAu"||. To this end, note that 
||W Tl || 2 = ||Vu" - Vv"|| 2 < 

Therefore, 

II Vu"|| < 2 + || Vv"|| < (^) 2 + M. 

We also have, for u™, the bounds 
fi 



f \\PAu n (;r)\\ 2 dr<C + C(t-t ), 
|K(-,r)|| 2 dr<C + C(t-t ), 



which are analogous to the bounds l|4.3|l and (|4.4|l for u, and can be proved by 
analogous arguments. Therefore, by Lemma 14. II 



supe"' / eT\\u?(;T)\\ 2 dT<C, 
t>0 Jo 

supe"* / e T ||PAu™(-,T)|| 2 dr < C. 
t>0 Jo 



Using these inequalities, one can show that 

-|||Vpu?f + S||Vu?f < C + C||PAu"|| 2 + C||u™|| 2 , (4.24) 

for all t £ [0,t*]. At this point, we need to restrict the time interval, since 

the constant C depends also on sup ||Vu"||, and we can assure this term to be 

t>o 

bounded, uniformly with respect to n, only in the interval [0, t*]. 
Using inequality H4.24fl . it follows that 

K(-,*)||<C. (4.25) 

Finally, inequality l|4.25|l allows one to prove 

||PAu n (-,t)|| < C, (4.26) 

for all t £ [0,t*]. We do not give the details of the proof, since it is completely 
analogous to the proof of Lemma f4. 41 □ 

Now we bound, for later use, the term VP n (ip n — £) = Vip n — VP„£. First, 
note that v™ satisfies 

(pn t , r ) + (Vv", V0") + Gem • Vu, 0") = (pf, </>"), (4.27) 
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for all <p n of the form 4> n (x) — dkW k (x). Subtracting equation l|3.8|) from 

fc=i 

equation l|4.27|l . we get 

(p«^,0") + (VV",V0 n ) = (p"e t n .0") + ((p-P n )(u t + u-Vu-f),0 n ) 

+ (p n u • Vu, 0") - (p"u" • Vu n , 4> n ). (4.28) 

On the other hand, taking the inner product of the fist equation in (|3.5|1 with 
<p n and integrating by parts, one gets 

(ft t , 4> n ) + (pu ■ v|, 4> n ) + (M ■ Vu > 4> n ) + (M • v|, <t> n ) + (V|, v</>") = 

((p-p n )(iK +u- Vu-f),0"). (4.29) 
Subtract equation l|4.29|l from equation (|4.28|l to get 

(p"0 t ,</>") + (V0,V0") = ((p-p")(u t + u-Vu-f),0") + ((p-p™)| t ,0 n ) 
+ (p"e£, 0") + (p"u • Vu, 0") - (p"u» • Vu", 0") 
+ (pu • V|, 0") + (p| ■ Vu, 0") + (p| ■ V|, 0"), (4.30) 

where 6 := ip n - |. Now, since P„0 = P„(i/>" - 1) = ^ n - P n |, one has 

(p"u-V0",0 n ) = (p"u-VP n 0,0") + (p"u-VP„|,<n, 
(pV*-Vu,0 n ) = (p^P„0-Vu,0") + (p"P n |-Vu,0"), 
(p"</>" • V^ n , 0") = (p>" • VP„0, 0") + (/>"</>" • VP„|, 0"). 

Therefore, 

(p"u • Vu, 0") - (p n u" • Vu", 4> n ) = (p"i/>" • Ve", <f> n ) + (p"e n • Vt/>", 0") 
+ (p"u • Ve", 0") + (p"e" • Vv", 0") - (p"u • VP„0, 0") - (p"u • VP„|, 0") 

- ( P n P n e ■ Vu, 0") - (p"P„| • Vu, 0") - (p>" • VP„0, 0") - (p"t/>" • VP n |, 0"). 

Moreover, since | = P„| + Q n £, one can show, after some computations, that 

(p"u • Vu, 0") - (p™u" ■ Vu", 0") + (pu ■ V|, 0") + (p| ■ Vu, 0") + (p| • V|, 0") = 
(p"t/>" ■ Ve", 0") + (p"e" • V0 n , 4> n ) + (p"u • Ve", 0") + (p"e" • Vv", 0") 

- (p"u • vp n e, 0") - (p n p n e ■ vu, 0") + ((p - p")u • vp„|, 0") 

+ ((p - p")P„| • Vu, 0") + (pu • VQ„|, 0") + (pQ„| • Vu, 0") - (p"t/>" • VP„0, 0") 

- (p"P„0 • VP„|, 0") + (p"P„| • VQ„|, </>") + (p"Q„| • VP„|, 4> n ) 
+ (p"Qnl • VQ„|, 0") + ((p - p")| • V|, 0"). 
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Applying this identity to 1)4. 30|) . and taking n = P n 9t, one obtains 



7"0t\\ 2 + \j t \\^P n 8\\ 2 = ((p-p")(u t + u-Vu-f+£ t + u-VP„£),P n t ) 

+ ((p - p n ){P n £ -Vu + f V£), PnOt) - (P n ^, Qn€t) (4-31) 

+ (p»e?, P„0 t ) + • Ve", P„0 t ) + (p"e n • V^ n , P n t ) 

+ (p n u • Ve", P n 6 t ) + ( P n e n ■ Vv", P n 6 t ) - (p n u ■ VP„0, P n 6 t ) 

- ( P n p n e • vu, p n o t ) + (pu • vq„£, p n 6> 4 ) + (pQ n $ • vu, p n e t ) 

- (p n *p n ■ vp„0, p n e t ) - ( P n p n e ■ vPni, P n e t ) 

+ {p n Pn£ ■ VQ„€, Pnflt) + 0"Qn£ ■ VP„£, P n t ) 
+ (p n Qn£-VQ n S,P n Ot). 

Now, we bound each term on the right hand side of identity l|4.31|) . Given e > 0, 
we bound 



\(o n 8t Qr,£ t )\ 


< 


c 


0/ 

"t 


2 + 


\(p n e? : P n O t )\ 


< 


e 


Ot 


2 + 


\{p n ip n -Ve n ,P n d t )\ 


< 


c 


0t 


l 2 + 


\{p n e n -Vip n ,P n e t )\ 


< 


c 


0* 


l 2 + 


\(p n u-Ve n ,P n 9 t )\ 


< 


e 


Ot 


2 + 


\(p n e n -Vv n ,P n O t )\ 


< 


c 


0* 


2 + 


\( P n u-vp n e,p n e t )\ 


< 


e 




l 2 + 


\{p n P n 0-Vu,P n t )\ 


< 


c 


o t 


l 2 + 


\(pu-vQ n ^p n e t )\ 


< 


c 




l 2 + 


\(pQ n $-Vu,P n d t )\ 


< 


e 


e t 


2 + 


\( P n xp n -wp n 9,p n e t )\ 


< 


e 


o t 


2 + 


\(p n p n e-wp n tP n e t )\ 


< 


e 


e t 


l 2 + 


\{p n P n £.-VQ n £.,P n 6 t )\ 


< 


e 


e t 


l 2 + 


\(p n Qn£-yPntP n e t )\ 


< 


e 


e t 


2 + 


\{p n Q n £.-VQ n £.,P n 6 t )\ 


< 


e 


e t 


l 2 + 



liv^f, 

l|Vu t || 2 , 

||PAt/,l 2 , 
||PA^"|| 2 , 



g(e) 
C(e) 

A n +1 

g(e) ) 

C(e)||VP„0|| 2 , 

•C(e)||VP„0|| 2 , 
C(e) 



An+l 

C(e) 



PA£|| 2 , 
ll^ll 2 , 



C(e)|!PA^"|| 2 ||VP„0|| 2 , 
C(e)||PA,£|| 2 ||VP„0|| 2 , 
C(e) 



g(e) | 

A n +1 

g(e) | 

An+l 



|PA£|| 4 , 
|PA£|| 4 . 
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It remains to bound |(7rg n , P n 9 t )\, where tt := p — p n and 

g™ := u { + u- Vu - f + £ t + u • VP„£ + P n £ ■ Vu + £ • V£. We begin by bounding 
g". 

Lemma 4.10. For all p, 2 < p < 6, the bound 

||g n M)HL < C + C||PA£(.,i)|| 2 + C||PA£(.,i)|| 4 

+ q|V« t (.,t)|| 2 + C||V^(.,t)|| 2 (4.32) 

holds for all t > 0. 

Proof. Since 2 < p < 6, we have 

l!g"M)HL < c{\\u t (;t)\\l P + \\u-vu(;t)\\h + WK-Mh + Hti-Mh 

+ ||u • VP„£(-,i)|| 2 P + \\P n £ ■ Vu(-,t)||| P + u • v€(.,t)||i,} 

< C{||Vu 4 (.,t)|| 2 + ||PAu(.,<)|| 4 + ||Vf(-,<)|| 2 + ||V^(-,<)|| 2 

+ ||PAu(.,i)|| 2 ||PA£(.,t)|| 2 + j|PA£(.,<)|| 4 } 

< C + C\\PA^M 2 +C\\PA^,t)\\ 4 + C\\Vn t (-,t)r + C\\VU,t)\\' 

□ 

Lemma 4.11. 7/6 < po < °°; fAen £/ie bound 



H;t)\\ir < C\\n{;t )\\lr + C(t-t ) ||VP„0(-,r)|| 2 dr 

+ T^-(t-io) 2 + ^-(i-io) f \\PAH(;r)\\ 2 dT (4.33) 

/io/rfs /or all t > and all r, 2 < r < — J/p n = oo, i/ien £/ie bound fZ! 551) 

6+Po 1 ' 

is vaZzd /or all t > and all r, 2 < r < 6. 

Proof. First note that 

7T, + u" • Vtt = (u™ - u) • Vjo. 

Since u = u + £, we write the equation above as 

TT t + u n ■ Vtt - (P n - Q n i - e") • Vp. (4.34) 

Let r belonging to the suitable interval depending on the value of pq. Multiply 
equation (|4.34(l by |7r| r_1 and integrate to get 

~IK,t)||£ P < / (P„0 - Q n £ - e") • Vp^y^dx 
r dt J Q 

< [Jj(P n O-QnZ-e n )-Vp\ r dxy U n W\ r dx 



< \\(Pne-Qn$-e n )-Vp\\^h\\ r L 



r-1 
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Therefore, 



■^\H;t)\\ L r < \\(P n e-QnS-e n )- Violin < \\P n - Q n £ - e n \\ LP \\Vp\\L?o, 

where p is chosen as — = if 6 < po < oo, and as p — r if po — oo. Note 

p r p 

that in the case 6 < po < oo, this choice of p implies 2 < < p < 6. In 

Po - 2 

the case po = °Oi we have 2 < p < 6. In both cases, p £ [2, 6] and we bound 



Jt h(;t)\\Lr < C(\\VP n e(;t)\\ + \\VQ n Z(;t)\\ + \\Ve n (;-t)\\) 



< c[\\vp n e(;t)\\ + 



c 



c 



-||PA^(.,t)|| , 



(A n +l) 2 (A n +i) 

where, for the last inequality, we used Lemma f2. II and inequality (|4.19|) . Inte- 
grating this inequality from to to t, 



H;t)\\lr < C{\\n(;t )\\lr + 



f \\VP n O(.,T)\\dT 
Jtn 



1 



to (A 



-dr 



n+1, 



la 



(AnH 



-||PAC(.,r)||dT 



< C^||7r(.,*o)||ir + (t-tb) / \\VP n 6(;T)\\ 2 dT 

J t 

+ J—(t-t ) [ ||PA£(-,T)|| 2 dr), 

A n+l Jto J 



•\n+l 



which is the desired bound. 

Getting back to inequality 14.31(1 . we have 



□ 



2 , / 1c ,„ fl „ a , C(e) ||V7 , „ a , C(e) 

An+l 



|V€ t | 



l|Vu*| 



2 , C{e) 



A n +1 A n +1 

PA€||* + C(e)||7r||i,||g" 



An+l An+l An+l 

+ C(e)||VP„0|| 2 {1 + \\PAi{, n \\ 2 + ||PA||| 2 } 



where, in the case 6 < po < oo, the inequality above holds for each r G 



3. 



6po 



6 +Po 



with p 6 



Po ~ 3 



,6 



chosen such that — I — = -. In the case 
r p 2 

1 1 1 



Po = oo, it holds for all r € [3, 6], with p G [3, 6] chosen such that 



Now fix e = — 
15 



1 

a 

2 



r p 2 

Integrating the inequality from to to t, and using 
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Lemma T4. 61 Lemma f4. 91 and inequality 1)4. 18|) . we get 

||VP„0M)|| 2 + f \\e t (;r)\\ 2 dr < \\VP n 8(;to)\\ 2 + -£-(t - to) + ^- f || V£ t (-, r)\\ 2 dr 

Jt A " + l A n+1 Jto 

+ f : ||Vu t (-,T)dr|| 2 + C nW(-,r)\\lr\\g n (-,r)\\ 2 LP dr + C f || VP„0(-, r)|| 2 dr. (4.35) 

A n +1 Jto Jto Jtn 



'to Jto Jto 

Adding inequalities (|4.33|1 and l|4.35|l . one obtains 



||VP„0(-,i)|| 2 + ||7r(-,t)||i r + f \\e t (;r)\\ 2 dr < || VP„0(-, t )\\ 2 + C\\n(; t )f Lr (4.36) 

Jt 

/ t ||V^(.,r)|| 2 dr + -^ r||Vu t (.,r)rfr|| 2 

Jtn A„_|_l J t 



A n +1 A n +1 A n+ i J tn A„+l J t(] 

+ C f ||7r(-,r)||i,||g n (-,r)||| P dr + C f ||V ■p n 0(-,r)|| 2 dr + C(t - t ) f \\VP n d(;r)fdr. 

J to J to J to 

Fixing i > to, and using Corollary 14 . 51 and Lemma ET71 one concludes that 



Mo)ll 



||vp„0(-,i)|| 2 + |k(-,t)||£,+ f \\e t (;T)fdT< \\vp n e 

Jt 

+ C\\n(-, t )\\h + + -^—{t -t ) + ^(t- t ) 2 (4.37) 

A n +1 An+1 A n +1 

+ C /* ||7T(-, r)|| 2 .||g"(, r)|| 2 P dr + C Al + (t - i O )}||VP n 0(-, r)|| 2 rfr, 

Jto Jto 

for all i S [t ,i]. Let A(t) := || VP„0(-, i)|| 2 + + / ||0 t (-, r)|| 2 dr. 



Thus, inequality 14.37fl gives 

A(t) < CA(t ) + -r— + T——(t - to) + T— (* - to)' 
A n +1 An+1 A„ + i 



+ C f {l+i-to + ||g n (-, r)|| 2 P } A(r)dr. 

./to 



Applying a corollary of Gronwall's Lemma(see [J, page 90, corollary 6.2), one 
gets 

A(i) < (cA(t ) + J— + ^—(t-to) + -^—(t- t ) 2 ) exp (c f (1 + t - t„ + r)|| 2 P )dr 
\ A n +i A„+i A„ + i / ^ J tQ 

We summarize the results in the following lemma. 

Lemma 4.12. Suppose ||V^"( , ,t)|| < K/Ak+i holds for a constant K > and 
all t in a given interval < to < t < t. Let £ as in problem and the 
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functions ir, 6, g™ defined as before. If 6 < po < oo then, for all t £ [to,t], one 
has 

\\X7P n O(;t)\\ 2 + h(;t)\\h + f \\0t(;r)fdT (4.38) 

J to 

< C (\\VP n 0(;t o )\\ 2 + ||7r(-,to)||ir + y~ !~ T^(* - *o) + T~~ — - M 2 ) e C/ 4 a(T)dT , 



/or r € 



6po 



6+Po 



, and p £ 

2 



3po 



,6 



chosen such that — | — = — , where 
r p 2 



a(t) := 1 + 1 — to + ||g™(-, 01 1 z,p ■ -tf Po — °°; then the bound holds for all r £ [3, 6] 

and p £ [3,6] sucft £/ia£ — I — = — . 

r p 2 

Remark 4.13. From inequalities (|4.32() . I|4.17|l . i|4.18[l . and Lemma f4. 61 one can 

bound 

/ a{r)dT < C + C(t - t ) + C(i - t )(t - t ). 

Therefore, 

exp | jf o(r)dr| < cxp {(5 + C*(i - to) + C(t - t )(t - *o)} • (4.39) 

From now on, we fix the constants C and C appearing in inequalities l|4.38[) 
and (|4.39() . We prove now that the bound for ||Vt/> n || required in lemma (|4.9|l 
and lemma (|4.12(l hold for n large enough. 

Proposition 4.14. There exist K > and N £ N such that if n > N , then 
\\Vip n (-,t)\\ 2 < for allt> 0. 

Proo/. Choose T such that M|(F(T)) 2 < -. Let X := 8C(1 +T + T 2 ) exp{<5 + 

~ ~ K 

CT + CT 2 } and let N to be large enough such that < S if n > N. Under 

A n +1 

these conditions, we have 



l|v^"(-,t)|| 2 < 



K 



(4.40) 



for all t > 0. Indeed, suppose that inequality 1)4. 4U[) does not hold. Thus, there 
exist n > N and t* > such that 



l|vv"(-,OH 2 = 



if 

An+1 



(4.41) 



We may have either t* < T or t* > T. If t* < T, consider t = 0, £ = 0, rj = 0, 
i — t* . In this case, ||VP„#| = ||Vi/>"j|. Therefore, using Lemma f4. 121 we have 



||V^ n (-,t*)r + ||7r(.,OII^+ MMVdr < 



c 



An+1 A„ + 

K K 

< 



c T C 



+ 1 



rp2 \ e C+CT+C'T 2 



8A„+i A„+i 
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which contradicts jOJl , If t* > T, apply Lemma Q2I with t = t*,t =t*-T 
and £(x,t), r)(x,t) satisfying 

T){x,t ) = p n (x,t ) - p{x,t ), 



to get 



||VV n (-,t*) - VP n £(-,t*)|| 2 + ||p(.,f) - + »,(■, t*)|| 



IK(-,r)-€t(-,r)|| 2 dr< 



c 



c 



c 



. j.2 \ e C+CT+CT 2 



Therefore, 

||VtT(.,OH 2 < 2(||v^ w (- ) **)-VP„€(-,**)ll 2 

K 



ivp n €(-,mi 3 ) 



< 2 

< 2 



8A, i+ i 
K 



Mt\\V£{;t )\\ F{T)' 



K 



3 K K 

< 



s 8A n +i 4A„+i / 4 A„+i A„ + i 
which again contradicts 1)4.4111 . This proves the proposition. 

5 Proof of Theorem EH 

Using the bounds (|4.19|) and 14.40fl . we directly get 

||Vu(.,t) - Vu"(.,t)|| 2 < 2(||V^(-,*)|| 2 + ||Ve"(.,<)|| 2 ) < 



□ 



C 



An 



+i 



(5.1) 



which is the first bound stated in Theorem 13.31 In order to prove the bound 
13.12f> for the density, first note that 



p t + u • Vp = 
p? + u" • Vp" = 0. 

Subtracting equation (|5.3|) from equation (|5.2I) . one gets 

(p - p n )t + u" • V(p - P n ) = (u« - u) • Vp. 



Now, if 6 < po < oo, let r £ 
111 



6po 
6 +Po 



Choose p £ 



2p 
Po-2 



,6 



(5.2) 
(5.3) 

(5.4) 
such that 



r p Po 
we obtain 



Multiplying equation l|5.4|) by \p — p n \ r 1 and integrating over fi, 



= / |p-p n | r_1 (u n -u) ■ Vpdx 

< ||p-p n ||27 1 ||(u w -u).Vp|Ur 

< ||p-p"||27 1 ||Vp|| L ™||u»-u|| iP , 
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If po = oo, the bounds above hold for all r € [2, 6] and p = r. Thus, 

±\\ p (.,t)-p n (.,t)\\ L r < G||u n (-,t)-u(.,t)|| L r < C\\Vu n (;t)-Vu(;t)\\. (5.5) 
Integrating inequality Ij5.5|l from to t and using 1)5. one gets 

\\p(-,t) - p n (;t)\\Lr < — ^-rt+ \\p(;Q) - p n (-,0)|| L , = -^-rrt, (5.6) 

(A„+l)2 (K + i) 2 

since p n (x,0) = po{x). This finishes the proof of Theorem 13.31 

A Proof of Lemma 14.11 

Suppose h(t) integrable, nonnegative and satisfying, for all t, t with < t < t, 

f h{r)dT < ai(t-to) + a 2 , (A.l) 

where a\ and a 2 are nonnegative constants. We first consider the case < t < 1. 
In this case, 



e" 4 



/ e T h(r)dT < f h(r)dT < f h(r)dT < a x + a 2 . 
Jt„ Jt„ Jo 



/to Jt JO 

Now, if t > 1, let n € N and r e [0, 1) such that t = n + r. Then, 

pt pn-\-r n pi pn+r 

I e T h(r)dT = / e T h{r)dT = ^ e T h{r)dT + / e T h{r)dT 
Jo Jo j* = i/j— i J n 

n pi pn+r 

< £V / h(T)dr + e n+r h{r)dT 

j—l Jj—l Jn 
n pj pn+1 

< £V / h{T)dT + e n+1 h{r)dT. 

j—l Jj—l Jn 



3 

Inequality l|A.l|) implies 



h{r)dr < a\ + a-2 , j = l,...,n + l. 

3-1 

Therefore, 

pt pn+r n+1 

e~* / e T h(r)dT = e - n ~ r / e T h(r)dT = eT n ~ r V e j / h(r)dT 
Jo Jo - =1 Jj-l 

< (oi + a 2 )e~ n - r Y e? = (oi + a 2 )e-"- r £ " +2 T - 
^— ' e — 1 

3=1 



(at + a 2 ) < (ai + a 2 )' 



e 2 



e - 1 ~ v y e - 1 
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Thus, 

supe * / e r h(r)dT < (oi + a 2 ) < oo, 

t>o Jo e - 1 

which proves the lemma. 
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